Two types of topological transitions in finite Majorana wires 
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Motivated by the recent advances in studying Majorana states in nanowires under conditions of 
superconducting proximity effect, we address the correspondence between the common topological 
transition in infinite system and the topological transition of other type that manifests itself in the 
positions of the poles of the scattering matrices. We establish a universal dependence of the pole 
positions in the vicinity of the common transition on the parameter controlling the transition, and 
discuss the manifestations of the pole transitions in the differential conductance. 



Majorana bound states have been predicted to exist 
in various condensed matter setups: 5/2 FQHE state 
[l[ , in vortices found in p + ip superconductors Q and 
in specific models of Id superconductors Q. The impor- 
tance of the Majorana states for quantum computation 
Q has brought them to the focus of the condensed mat- 
ter research [5j. Next step were the suggestions to realize 
the Majorana states in more experimentally feasible se- 
tups, those include topological insulators [1,0] and semi- 
conductor nanostructures with big spin-orbit interaction 
brought in proximity to s-wave superconductors. Two- 
[H, Q and one-dimensional 13, U| 
been considered. 



nanostructures have 



The observation of Majorana bound states in Id 
nanowires has been reported by several groups by mea- 
suring zero-bias conductance peak fl2l - [l4l ] and Air Joseph- 
son effect [lj|. The signature of Majorana's is their emer- 
gence upon a topological transition [l6| separating the 
regions of parameter space with and without zero-energy 
excitations. In all cases the experiments have been per- 
formed with finite and rather short wires. This brings 
about the question: how a topological transition taking 
place in infinite system is manifested in properties of a 
finite wire. 

Strictly speaking, this common topological transition 
is absent in a finite system where excitation energies con- 
tinuously depend on the control parameter of the tran- 
sition and are never precisely zero [ 3| . llCl Il7 |. This 
may be shown in several ways. In [17| we gave the 
most general formulation in terms of the topology of the 
energy-dependent scattering matrix characterising a fi- 
nite nanostructure. Same study revealed a topological 
transition of other kind that takes place in finite systems 
and manifests itself in the properties of the poles of the 
scattering matrix. The topological number in this case is 
the number of poles at purely imaginary energy, and the 
topological transition is the change of this integer even 
number upon the continuous variation of the control pa- 
rameter. 

In the present work we link these two topological tran- 
sitions of different types: bulk one and finite system 



one. We show that in general case the common topo- 
logical transition is accompanied by the pole topological 
transition[l8j. The points of the transitions differ at the 
scale inversely proportional to the wire length. We im- 
plement the generic model of the Majorana wire that is 
always valid in the vicinity of the transition point and ob- 
tain the universal dependence of the pole positions on the 
control parameter and a single parameter characterizing 
the coupling of the wire to a normal metal lead. We dis- 
cuss how the same correspondence occurs for more spe- 
cific models and how the universal picture is manifested 
in a transport measurement. 

The "standard" model describing a Majorana wire en- 
compasses a single-band spectrum that includes spin- 
orbit interaction, proximity effect from the bulk super- 
conductor and spin magnetic field [ToL 11 1. Let us derive a 
phenomenological effective model valid near the common 
topological transition point. We can start with a multi- 
mode wire where the spectrum at each k is described by 
a general Hamiltonian matrix H(k) in the space of the 
modes and Nambu index. The general symmetry of BdG 
equations requires H(k) = —H T {~k) in a certain 
(Majorana) basis. The common topological transition 
takes place when an eigenvalue of H(k = 0) passes in- 
dicating a closing of the superconducting proximity gap 
in the wire. 

Owing to BdG symmetry, the zero eigenvalue is doubly 
degenerate. Thus we concentrate at the two modes cor- 
responding to the eigenvalue. Near the transition point 
the general form for this Hamiltonian in Majorana basis 
reads H(k = 0) = acr y , where u's here and below are 
usual Pauli matrices. The phenomenological parameter 
a controls the proximity to the transition and is a func- 
tion of physical control parameters like magnetic field or 
chemical potential, a = in the transition point. Ex- 
panding near k = and taking into account the BdG 
symmetry, we find two possible terms oc ka x and cx ka z . 
The combination of the two can be brought to oc ka z by 
a rotation of the pseudospin about y axis. This brings us 
to the generic Hamiltonian we will use in further consid- 
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cration: 



H = vko~ z 



(1) 
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It has been first introduced in 

Let us turn to a finite wire The boundary conditions at 
the wire ends must be consistent with the current conser- 
vation. The operator of current reads I = v<r z so the con- 
servation implies that the wavefunction "J = {ipi,ip2} T 
has to satisfy: 

l^il 2 = h/> 2 | 2 , (2) 

Near zero energy the wavefunction is real, and we are 
left with binary choice ip\ = ±-02 ■ We fix the signs to + 
at the right end of the system and — at the left one. In 
this case, in the limit of infinite wire length the Majorana 
states are formed at a < while the phase at a > is 
topologically trivial (Fig. QJ>). 

Let us now contact the left end of the wire with a nor- 
mal metal lead and describe the situation in therms of the 
scattering matrix from/to normal lead modes. Scattering 
matrices are very useful objects to study the properties of 
the superconducting junctions [22T]. They incorporate rel- 
evant details of the setup in few parameters and allow to 
compute different properties of the junction, like conduc- 
tivity. They also allow for the topological classification 
of the junction in a concise way [17] . 

There are two interesting modes in the wire propagat- 
ing in opposite directions. The scattering matrix of the 
contact S c is in the basis of the incoming waves in the 
lead and the single mode of the wire and is thus of size 
M+lxM+1, M is the number of modes in the normal 
lead. We separate it into blocks as: 



rn r 2 i 
fi2 -r 



(3) 



Here fu is M x M matrix of the (Andreev)reflection to 
the leads that also incorporates the scattering in all other 
wire modes, r^i, 12 are scattering amplitudes from/to the 
wire and r is a number, which gives the reflection am- 
plitude in the wire (r = 1 corresponding to the wire 
isolation) [23| . By virtue of BdG condition r is real at zero 
energy. Since the interesting energy dependence comes 
from the wire, we can neglect the energy dependence of 
S c . 

To get the full M x M scattering matrix in the space of 
normal lead modes, we need to combine the S c with the 
scattering amplitude S w that describes the propagation 
along the wire, reflection from the right and the propa- 
gation back to the left end. This amplitude is easy to 
find from the Hamiltonian (1) and reads 



S-nr — e ^ — 



)th(Va 2 - e 2 L/v) + -0 



coth(Va 2 — e 2 L/v) 



(4) 
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FIG. 1: (a) The setup: the Majorana wire of the length L at 
the top of a superconductor is connected to a normal-metal 
lead. The total scattering matrix at low energy incorporates 
that of the contact (S c ) and energy-dependent scattering ma- 
trix describing propagation in the wire, S w . (b) Sketch of the 
(continuous) spectrum in the limit of infinite L: a Majorana 
level emerges upon the common topological transition, (c) 
The common topological transition becomes a crossover for 
finite L. The quantized energy levels (real parts of the pole 
energy positions) are sketched versus the control parameter a. 
The lowest level reaches at the point of the pole topological 
transition. Dashed lines give the imaginary parts of the pole 
positions for the lowest level. 



L being the wire length. The whole peculiarity of the 
limit of infinite wire may be seen from this formula. If 
we formally set L — > = sign(a). The ma- 

trix is thus energy-independent and topologically trivial 
(nontrivial) for a < (a > 0). 

The full scattering matrix thus reads 



S = fn +f 2 ie 



;ri2- 



(5) 



We concentrate on poles of this matrix those are solutions 
of 



\/a 2 — e 2 coth(\/ a 2 — e 2 L/v) + a — ie 



1 



1 + r 



0. (6) 



At finite length L the common topological transition 
becomes a crossover taking place in an interval of a of 
the order of effective level spacing in the wire v/L and at 
the corresponding energy scale. We aim to describe this 
universal crossover. To this end we rescale a, e to dimen- 
sionless units a = ay-, e = ey. The equation becomes 



Va 2 - P cotl^v/^ 2 - ~e 2 ) + a - "3— j; = °- ( 7 ) 

Numerical solutions for pole positions are shown in Fig. 
EH,b for two values of r as functions of the control param- 
eter a. We see a sharp feature in the crossover region: 
the pole topological transition. At this point, the real 
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part of the energy of the lowest pole becomes strictly 
zero. This occurs at finite negative values of a. The 
higher the transmission through the S c , the closer to 
is the transition point. This dependence is presented in 
FigHt- In the limit of low transmissions, the pole transi- 
tion takes place at \a\ ~ ln(l — r) where the exponentially 
small splitting of Majorana states matches small decay 
rate of the left-end state to the normal metal. The real 
parts of energies of all other poles follow the hyperbola- 
like curves indicating formation of discrete energy levels 
in the wire above the gap edge \a\. The same transi- 



tion is seen in imaginary parts of energy positions as a 
splitting of the curve corresponding to the lowest pole. 
The upper (lower) parts of the split curve give the decay 
rates of the left(right) end Majorana state. The decay 
rate for the Majorana "buried" at the right end falls off 
exponentially with increasing |a|: e ~ 2iaexp(— 2|a|) jzS- 
Let us address the dependence of the pole positions on 
a near the pole transition point. For this, we expand 
near the transition point a = ao, e = ieo to obtain the 
relation between the deviations 5t,5a from the transition 
point in the lowest non- vanishing order: 



-ag(l + 2e Q ^) - 2elii (-1 + n 2 ) + 2a e 2 (l + gg (1 - 3^ 2 )) + qgep (-3// - £o (-1 + 4e ^ + n 2 )) 
2e (a 2 + el) (a (l + 2a Q )n + e (-1 + a (-1 + m 2 ))) 
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FIG. 2: The real (a) and and imaginary (b) parts of the pole 
energy positions versus the control parameter a at two values 
of the reflection amplitude r: r — 0.86 (green dots, almost 
isolated)and r — 0.34 (red crosses, almost transparent). The 
part of (a) within the rectangular is replotted in (c). The pole 
topological transition occurs at ao = —1.2 for r = 0.86 and 
ao = —3.3 for r = 0.34. (d) The dependence of the transition 
point a on r. 



Here /i = jt^- This gives square root splitting of either 
real parts of the energy positions 5e = ±^/Sa/C at 8a > 
or imaginary ones, Se = ±iy/\5a\/C at 5a < 0. This 
square root dependence of 5e on 5a is in full agreement 
with Fig. [2k,b. 

The experimentally observable quantity is the differen- 
tial conductance of the contact, G(e — eV), V being the 
voltage drop at the contact. In terms of the scattering 
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FIG. 3: The universal bias-dependent contribution to dif- 
ferential Andreev conductance of the contact versus en- 
ergy/voltage. Solid, dashed, dotted curves correspond to 
r = 0.9, 0.5, 0.1 respectively, (a) 3 = 2, long before the tran- 
sition, (b) 5 = 0(c) 5 = —2 in the crossover interval, (d) 
a — —4 long after the transition. The Majorana at the far 
end of the wire is manifested as a narrow dip at zero bias. 



matrix, the conductance reads G = ^jTr(l — a y Sa y S^). 
Substituting S in the form of ((5|), we obtain a universal 
energy dependence of the conductance in the crossover 
region, 



G(e) = G + G 1 f(a,e); (9) 

(1-7-2)2 

f(a,e)= { > (10) 

1 + r jS + 2r cosx(a, e) 



The dependence is governed by the universal function 
/(a,e) (0 < / < 4) while the non-universal coefficients 
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Go, Gi depend on the details of the S c 

„2 



Go 



2irh 



Tr (1 - (r u 



.-1.S 



r 2 iri2jo- J/ r 11 <7j /; ) 



Gi 



e 2 Tr {cryf2if 1 2<Jy(f[ 2 f{ 2 + (r - 



)) 



2vrft 



(1 - r 2 ) 2 



(11) 



The coefficient Gi ~ e 2 /fi and can be of any sign while 
Go can be much bigger than e 2 /h. The function / (Fig. 
[3]) at any r gives a sequence of peaks associated with the 
poles of the scattering matrix. The peaks are narrow in 
the isolation limit r — >• 1. Before the transition, the peaks 
are far from zero energy. Upon the crossover, the peaks 
come close to zero and almost merge near the transition 
point. However, they never merge to a single peak: the 
Majorana state at far end of the wire is manifested in the 
conductance as a dip that becomes increasingly narrow 
upon increasing —5. 

Since the poles always have a finite imaginary part, 
and the conductance is defined at real energy, there is 
no singularity in /(e) at the point of the pole topological 
transition. However, this singularity can be in princi- 
ple identified from the experimental data by numerical 
analytical continuation to complex energy plane. The 
pole topological transition in principle does not require 
the common topological transition and can occur, for in- 
stance, in strongly disordeded wires 24J. The signatures 
of pole transitions have been investigated for the disor- 
dered wire model (l9j in the parameter region where no 
common transition can guarantee the presence of Majo- 
rana mode. 

Another setup proposed 
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of Majorana fcrmions encompasses normal leads at both 
ends of a finite nanowire. Also in this case the common 
topological transition is accompanied by a pole transition 
and proceeds in a similar way. The qualitative difference 
is that far below the transition both Majorana states re- 
tain a finite width and each of the two associated poles 
is manifested only in the scattering from the correspond- 
ing end of the wire. In the model under consideration, 
the Majorana splitting retains the same sign. More de- 
tailed models, e.g.[lo|. predict spectacular oscillations of 
the splitting [25[ ■ We stress that in the limit of the long 
wires Lkp 3> 1 such oscillations can only start far from 
the common topological transition, that is, at the values 
of the control parameter that are parametrically bigger 
than v/L. 

To conclude, we have formulated and studied a univer- 
sal model that describes the crossover in the vicinity of 
the common topological transition for finite clean Majo- 
rana wires. Importantly, we have shown that the sharp 
pole topological transition takes place in the crossover 
interval of the control parameter and computed the de- 
pendence of the pole positions on the contol parameter 
in this interval. We have also found a universal shape of 
differential conductance for this model, this enables its 
starightforward experimental verification. 



We stress the universal character of our conclusions, 
in particular, the predictions for the conductance: those 
should hold in any sufficiently long wire with small dis- 
order in the vicinity of the topological transition. Some 
features of our results have been seen in Ref. [12|: the 
authors have observed a narrow zero-bias peak on the 
background of a wider dip as seen in Fig. [3ji (assuming 
Gi is negative). From the other hand, no regular pat- 
ten of peaks moving to zero upon changing the control 
parameter has been observed so far. More experimental 
data, in particular, for longer wires are required to clar- 
ify the discrepancy that can be due to sufficiently strong 
disorder or finite temperature effects. 
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